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ABSTRACT 

We present an analysis of BeppoSAX observations of the unique transient bursting X-ray pul¬ 
sar GRO J1744-28. The observations took place in March 1997 during the decay phase of the 
outburst. We find that the persistent broadband X-ray continuum of the source is consistent 
with a cutoff power law typical for the accreting pulsars. We also detect the fluorescence iron 
line at 6.7 keV and an absorption feature at ~ 4.5 keV, which we interpret as a cyclotron line. 
The corresponding magnetic field strength in the line forming region is ~ 3.7 x 10 11 G. Nei¬ 
ther line is detected in the spectra of the bursts. However, additional soft thermal component 
with kT ~ 2 keV was required to describe the burst spectrum. We briefly discuss the nature of 
this component and argue that among other possibilities it might be connected with thermonu¬ 
clear flashes at the neutron star surface which accompany the accretion-powered bursts in the 
source. 
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1 INTRODUCTION 


The peculiar transient X-ray source GRO J1744-28 also known 
as the bursting pulsar, was discovered on December 2, 1995, with 
the Burst and Transient Source Experiment ( BATSE ), on board 
the Compton Gamma Ray O bservatory ( CGRO ) dFishman et al.l 
1 19951 : iKouveliotou et al.| [l996l). Three major outbursts have been 
observed since the discovery. The first lasted from the end of 
1995 until the beginning of 1996, the second from Decem¬ 
ber 1 996 until April 1997 , and the third from January to May 
2014 dDeeenaar et al.ll20l4 lYbunes et alJkOli iD'Ai et al .1120151) . 
The peak o utburst lumin osity in X-rays reaches L x ~ 10 37 - 
10 38 erg s -1 JWoods et ajj|l999| ), while the quiescent luminosity is 
~ L x ~ 10 33 erg s -1 dWiinands & Wandl2002l : iDaigne et~afll2002l : 
iDegenaar et al.l2012h . The source exhibits coherent pulsations with 
a period of ~ 0.467 s, associated with the rotation of t he ne u tron 
star i.e. is an accreting X-ray pulsar in a binary system. [Finger et all 
( 1996 ) determined the orbital period and semimajor axis to be 
P orb = 11.8 days and a = 1.12 R© (~ 7.8 x 10 10 cm) respectively. 
The mass function was estimated to be f(M) = 1.36 x 10 -4 M© 
implying, for the canonical neutron star mass of M NS = 1.4 M©, 
accretion via Roche lobe overflow from a strong ly evolved red gi ¬ 


ant remnant with mass about 0.2 - 0.7 M© (Daumerie et alj 


Miller 19961: Is turner & Dermerl 1 19961 : iBildsten & Brownl 


1996 


1997 


Rappaport & Jossl 1 19971) . The photoelectric absorption does not 


change significantly during the outbursts and is comparable with 
interstellar absorbtion N H ~ (5 - 6) x 10 22 cm -2 in the direction 
of the Galactic cen ter where the source is thought to be loca ted at 
distance of ~8 kpc dDotani et al.ll 19961 : INishiuchi et al.ll 19991) . 

Probably the most unusual property of GRO J1744-28 is 


the relatively short ~ 10 s bursts observed from the source. Since 
the discovery, several thousan ds of bursts have been detected, all 
at high luminosity dNishiuchi et al . 1999 ). During the bursts the 
flux increases by an order of magnitude, and it is often followed 
by a drop to below the pre-burst level for several tens of sec ¬ 
onds to minutes depending on the burst fl uence dLewin et al .11 19961 : 
INishiuchi et al.|[l999l : lYounes et al.ll2015l) . Pulsations are observed 
during the bursts, althoug h with a phase shift of ~ 5% with respect 
to the persistent emission dKouveliotou et al.fl996l : IStrickman et al.l 
fl996|) . 

A similar bursting activity has been observed only in another 
source the transient LM XB MXB 1730-335, usually referred to as 
the Rapid Burster (RB) dLewin et al]ll976l : lLewin et al.lll993l) . For 
both sources, it has been suggested that the origin of the bursts 
could be due to accretion flow instabilities which intermittently 
enchance the accretion rate onto the neutron star (so-called Type 
II bursts). However details are still unclear. For instance, unlike 
the RB, Type I bursts, associated with thermonuclear flashes on 
the surface of neutron star, are generally thought not to occur in 
GRO J1744-28. On the other hand, contrary to the RB case, the 
duration of the GRO J1744-28 bursts is remarkably stable (~ 10 s), 
and a large fraction of the bursts exhibit a temporal profile char¬ 
acterized by a fast rise followed by an exponential decay typical of 
classical bursters. Type II bursts in RB have rather irregular profiles. 
These differences are likely related to the different magnetic field 
strengths of the neutron star companion, which is expected, given 
the presence of pulsatio ns, to be stronger in GRO J1744-28. In the 
case of GRO J1744-28. lFinger et al.l d l996l) reported an upper limit 
on the magnetic field strength of B ^ 6 x 10 11 G. This estimate is 
based on the requirement that the plasma at the inner edge of the ac- 
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cretion disk (disrupted by the magnetosphere) moves faster than the 
magnetic field lines as otherwise the accretion would be ce ntrifu- 
gally inhibite d (so - called “propeller” regime. Iniarionov & Sunvaevl 
1 19751) . Later ICuil dl997h found a possible transition to the “pro¬ 
peller” regime in RXTE data at lower luminosity, which allowed 
for a n estimate o f the the magnetic field at ~ 2.4 x 10 11 G. Recently 
IP’Ai et aD 12015 ) reported the observation of an absorption feature 
at E~ 4.7 keV, during the 2014 outburst of the source. Interpreting 
the line as due to cyclotron resonance scattering, the authors esti¬ 
mated a magnetic field of ~ 5.3 x 10 n Gin good agreement with 
earlier estimates. 

In this paper we report on independent discovery of the same 
feature in BeppoSAX observations carried out during the outburst 
in 1997 (Sect. 3. 2fl We analyze three BeppoSAX observations of 
GRO J1744-28 carried out during the outburst in 1996-1997 and 
present a detailed spectral and timing analysis of the persistent flux 
and the bursts. We also briefly discuss the origin of the bursts in 
GRO J1744-28. The instruments and methods of the analysis are 
described in Section 2. Results from the timing and spectral anal¬ 
ysis are summarised in Section 3. In Section 4, we discuss the ob¬ 
servational results and finally provide a summary of the paper in 
Section 5. 


2 OBSERVATIONS 


BeppoSAX was an Italian X-ray astronomy satel lite, with Dut ch par - 
ticipation, which operated from 1996 to 2002 ( iBoella et~aPll997 ). 
It covered a broad energy band 0.1-300 keV with four Narrow Fiel d 
Instruments (NFIs) and two Wide Field Cameras (llager et al.l 19971) . 
The NFIs included the Low-Energy C oncentrator Spectrometer 
(LECS, 0.1-10 keV. IParmar et al.|[l997h . three identi cal Medium- 
Energ y Concentrator Spectrometers (MECS, 2-10 keV. lBoella et al.l 
Il997 ). and two collimated high energy detectors with good en¬ 
ergy resolution and low instrumental background: the High Pres¬ 
sure Gas Scintillation Proportional Counter (HPGSPC, 4-120 keV, 
FWHM e ne rgy re solution of 8% at 10 keV and 5.5% at 20 keV, 
iManzo et"ai1ll997l) and the Phoswich Detection System (PDS, 15- 
300 ke V, FWHM energy res olution of 24% at 20 keV, and 14% at 
60 keV, iFrontera et all 19971) . 

GRO J1744-28 was observed by BeppoSAX several times 
from March 1997 to April 1998 (Table [lj. In Figure Q] the long 
term source light curve observed by the All Sky Monitor onboard 
the Rossi X-ray Timing Explorer (ASM RXTE . iBradt et al.l 19931 
with superimposed BeppoSAX observations is presented. The first 
three observations of BeppoSAX were performed in the declining 
edge of the 1997 outburst, with a total exposure time for MECS 
of approximately 270 ks. In our work we use LECS, MECS, and 
PDS in 0.7 - 4 keV, 2-10 keV and 15 - 120 keV energy ranges 
respectively. The HPGSPC was, unfortunately, not operating dur¬ 
ing the first three observations. In the last three observations, the 
source was in quiescence and below the sensitivity thresholds of 
BeppoSAX , so it was not detected. The data from the quiescence, 
however, proved very useful for the background subtraction of the 
outburst PDS data. 

The standard BeppoSAX pipeline was used for the processing 
of the data. HPGSPC and PDS were operated in rocking mode to 
monitor local background and its variation along the orbit. The 


1 We learned of the content by D’ai et al when the manuscript was in the 
final steps of preparation. 


Table 1 . Observations of the X-ray source GRO J1744-28 by BeppoSAX. 
Luminosity is estimated using unabsorbed flux in the energy range 2-10 
keV (MECS) for the distance D = 8 kpc. 


MJD 

Exposure 
time, ks 

Period, s 

L x , 10 37 erg s 1 
2-10 keV 

Orbital 

phase a 

50528 

117.5 

0.467044(1) 

2.8 

0.9 

50534 

101.9 

0.467044(1) 

1.8 

0.4 

50550 

51.6 

0.46705(1) 

0.46 

0.8 

50905 

62.5 

source in 

0.8 

50911 

34.5 

quiescent state 

0.3 

50913 

67.5 



0.5 


a Ephemeris from lKouveliotou etaO ( 1996) 



MJD 


Figure 1 . The light curve of GRO J1744-28 observed by the ASM onboard 
RXTE. Dashed lines mark the BeppoSAX observations (Table[I}. 


source counts for LECS and MECS were extracted from source 
centered circles with radii of 4 arcmin, while the background from 
an annulus with outer radius of 18 arcmin. A detailed description of 
BeppoSAX data reduction procedures can be found in the handbook 
for BeppoSAX. 


3 DATA ANALYSIS AND RESULTS. 

3.1 Timing analysis. 

The 2-10 keV MECS light curve of the first observation (MJD 
50528) is presented in Fig[2] Bursts are clearly visible in the first 
two observations, whereas none is detected in the last observation 
carried out close to the end of the outburst at significantly lower lu¬ 
minosity. We detected 25 and 17 major bursts in the first two obser¬ 
vations respectively, with a flux increase of a factor of 10 or more. 
Typical burst duration and the interval between the bursts are ~10 s 
and ~ 2000 s respectively (see also Fig. [3] where all detected bursts 
are reported). The inter-burst to burst fluence ratio in the energy 
range 2 -20 keV is between 3 and 18, i.e. com patible with earlier 
reports (iLewin et al.lll996l : Ijahoda et"aPll997h . The burst profiles 
exhibit variety of shapes: most are rather symmetric and in this 
respe ct are similar to Type II bursts observed in RB (iLewin et al.l 
Il976h . However, about twenty percent (for example, the next to last 
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one in Fig. [3} exhibit the characteristic exponential decay typical 
of thermonuclear Type I bursts. The source also exhibited a few 
“smaller” bursts, a factor of ~ 6 less l uminou s than the majo r bursts. 
These have already been reported bv lNishiuchi et al. ( 1999 ) based 
on the ASCA data. We confirm the presence of the low luminosity 
bursts, although the statistical quality of the data is not sufficient 
for any detailed analysis. 

As reported earlier by several authors (iLewin et al ]&; 
lYounes et al .1120151) . a flux drop following the burst is a relative 
common feature in GRO J1744-28. Such a drop is observed in 
more than half of the bursts in our sample. In Fig. [2] a fragment 
of the light curve in the 2-10 keV is presented. The decrease of 
the mean flux from the source after the burst, for about 100 sec, is 
evident. Note that a similar flux depression is observed after Type 
II bursts of the RB and is usually associated with the depletion of 
the inner edge of the accretio n disk after the bursting event due to 
enhanced accretion ( iLewin et"al .119961) . 

To study coherent pulsations, we corrected the photon ar¬ 
rival time for Doppler delays due to the orbital motion of 
the spacecraft an d the pulsar (using the orbital ephemeris of 
iKouveliotou et al.l 19961). Using the phase -connection technique 
[see for instance!( Doroshenko et al.l [20 id) we measure the spin 
period to be P = 0.467044(1) s, P = 0.467050(1) s and P = 
0.46705(1) for the March and the two April 1997 observations 
(all uncertainties are for 90% confidence level unless stated oth¬ 
erwise). Folding the light curves reveals sinusoidal, single-peaked 
pulse profiles (Fig. 0]). We find that the phase of the pulsed pro¬ 
files during the bursts is shifted by ~ 5% relativ e to the per sistent 
state, which is consistent with earlier reports by IKouveliotou et al.1 
( 1996 ); lstrickman et al.l i 1996 b. The morphology of the pulse pro¬ 
files appears to be consistent between the two observations and 
shows no apparent dependence on luminosity. 

The fraction of pulsed emission defined as the ratio (F max - 
F m m)/(Fmax + ^min), where F max and F min are the maximum and 
minimum source flux does change with energy. To estimate it we 
assumed standard backgrounds for all instruments and also took 
the contamination from the nearby sources to the PDS data (see 
Section lT2l) . We found that the pulsed fraction increases between 
2 and 10 keV from about 5% to 20% (with a notab le exceptio n of a 
region close to the fluorescent iron line, see also lD’Ai et al.ll2015b . 
and then remains fairly constant between ~ 10% and ~ 30% for 
the first two observations and between 20 - 80% for the third one. 
From 0.1 to 2 keV, the pulse fraction decreases from 20 - 30% to 
5% although our findings are at the lowest energies and for the third 
observation are hampered by low statistics (Fig. [5]). 

3.2 Spectral analysis. 

GRO J1744-28 is located in the crowded Galactic Center region 
and a number of sources fall into the field of view of the tele¬ 
scope. This might potentially affect the spectral analysis, particu¬ 
larly for non-imaging instruments. Even the LECS and MECS ob¬ 
servations could be affected as the deep Chandra observation of 
the field reveals ~ 40 sources within the standard 4' extraction re¬ 
gion. However, GRO J1744-28 is by far the brightest among these 
sources even in quiescence, so for the imaging instruments MECS 
and LECS source confusion is not an issue. We followed, therefore, 
standard analysis procedure. 

The situation for the non-imaging instruments is more com¬ 
plicated. The MECS image reveals at least two sources besides 
GRO J1744-28 within the field of view of collimating instru¬ 
ments. These are detected irrespective of the brightness of the pul- 



0 5 10 15 20 25 30 35 


Time, hours from MJD 50528.293 

Figure 2. MECS light curve with bursts of the pulsar GRO J1744-28 in 
1997-03-21 with 2 s bin time in 2 - 10 keV energy range. On the zoom 
picture the flux drop after the burst is shown. 



Phase 


Figure 4. Pulse profiles of the pulsar GRO J1744-28 in 1997-03-21 for 
LECS (0.1-4 keV), MECS (4-7-11.5 keV) and PDS (12 - 45 - 100 
keV). The spin period is P = 0.467044(1) s. 
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0 10 20 0 10 20 0 10 20 0 10 20 0 10 20 0 10 20 

Time from the flare start, s 



0 10 20 0 10 20 0 10 20 0 10 20 0 10 20 0 10 20 

Time from the flare start, s 


Figure 3. Light curves of all burst of the pulsar GRO J1744-28 in 1997-03-21 (top) and 1997-03-27 (bottom), MECS. 


sar (see Fig.[6]> and are plausibly identified as IE 1743.1-2843 and 
Sgr A*. To account for the in-orbit instrumental background which 
is rapidly changing and the for the contribution of the sources con¬ 
taminating the PDS data, we extracted the background-subtracted 
spectra for all observations (including those in quiescence) using 


the standard pipeline. The standard pipeline allows to account for 
instrumental in-orbit background. Then we subtracted quiescent 
spectra (obtained during the 1998 observations) from the source 
spectra obtained for the first three observations. Apparently, an 
important assumption here is that the flux of the contaminating 
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Figure 5. Pulse fraction of the pulsar GRO J1744-28 for 1997-03-21 - solid 
black line, 1997-03-27 - dashed red line, 1997-04-12 - dashed dotted blue 
line. 




Figure 6. MECS image centered at GRO J1744-28 in quiescence on 1998- 
04-02 overlaid with contours from an observation in bright state (taken on 
1997-03-21). Two additional sources detected in the MECS field of view 
contaminate also the PDS spectrum. 


sources remains constant. To verify that we compared the soft X- 
ray fluxes measured by the MECS which indeed remained constant. 

F or spectr al fitting, we used the XSPEC package (version 
11.3.2. lArnaudll 1996h . We found that the observed broadband con¬ 
tinuum spectrum can be described with several phenomenological 
models generally used for the spectra of accreting X-ray pulsars: 

1) a power law plus a high-energy cutoff model, where the 
transition is smoothed with a Gaussian 

PHC F~ r I 1 (E ^ Ecut) m 

y e -(E-E cut )/E Md {E > ^ ’ 



E, keV 


Figure 7. The best-fit residuals for spectra derived from the observation 
1997-03-21 modeled as described in the main text. The best residuals ob¬ 
tained with the power law plus high-energy cutoff model, PHC(E), at the 
top panel. 


2) a power law with Fermi-Dirac cutoff: 

FDCO ~ E~ r (1 + e ( E ~ E **)/ E *M J- 1 ; 

3) the Negative and Positive po wer law Exponential model 
dMiharalll995l : lMakishima et HEwl): 

NPEX = (AiE~ a ' +A 2 E +a2 )e~ E,EfM ; 

4) the CompTT model by Titarchuk dl994h describing the 
comptonization of soft photons with temperature kT s in a hot elec¬ 
tron plasma with kT e and a hot electron optical depth r e . 

In all cases the inclusion of the fluorescent iron K a line at 
~ 6.7 keV (modelled with a simple Gaussian profile) was also re¬ 
quired. All continuum models give a reasonably good x^ ed (from 
0.94 for the PHC model to 1.25 for CompTT model, see Table [2}. 
To demonstrate the evolution of the parameters with the luminosity 
we used PHC based on the formal statistical quality and stability 
of the fit (Table [3}. Regardless of the assumed continuum model, 
some residuals below 1 keV, and around 4-5 keV (Fig. ID were 
clearly observed. The absorption-like feature at ~ 4 keV was mod¬ 
elled with a multiplicative line with either a gaussian or lorentzian 
profile. This significantly improved the fit (see residuals of Fig. [8}. 
The parameters of the line for the gaussian profile are presented in 
Table |2] For the lorentzian profile (CYCLABS in XSPEC) we get 
comparable values, i.e. E cyc = 4.0 + 0.2 keV, cr = 1.3 + 0.2 keV. 
The chance probability for fit i mprovement based on the f-test is 
marginal ^ 10 -30 (see, however. IProtassov et alJ2003) . The feature 
also appeared in the phase resolved spectra (see Fig. [9}. Our find¬ 
ings nicely agrees wi th the absorption-like feature recently reported 
bv D’Ai et al.l d2015h based on XMM-Newton and INTEGRAL ob¬ 
servations of the 2014 outburst. We note that lYounes et al.l d20 1 5h 
detected an absorption feature at lOkeV which is not observed in 
our data, although the lack of HPGSPC data makes it difficult to 
draw a definitive conclusion. 

We have also performed a pulse phase resolved spectral anal¬ 
ysis to study the variation of spectral parameters with the angle of 
view to the neutron star. Only data from the observation on March 
21, 1997 with highest luminosity and best counting statistics was 
considered. To describe the phase resolved spectra, we used the 
same model as for the phase averaged spectrum, i. e. a power law 
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1 10 100 
E, keV 


Figure 8. The best-fit unfolded average persistent spectrum, modelled with 
the PHC continuum is shown in the top panel. We also present the best fit 
residuals without and with the inclusion of a cyclotron line at ~ 4.5 keV. 
The best fit spectrum of the bursts modelled with PHC and a black body 
continuum is also presented together with the residuals for the observation 
1997-03-21. 


Table 3. Spectral parameters of the persistent spectrum low mass X-ray 
pulsar GRO J1744-28 observed by BeppoSAX modeling by the PHC with 
smoothing gaussian continuum model. All the energies and the line widths 
are given in keV. 


Parameter 

1997-03-21 

1997-03-27 

1997-04-12 

N a 

5.7(2) 

5.7(2) 

5.4(1) 

r 

1.26(7) 

1.29(5) 

1.24(3) 

Fcut 

18(1) 

19.3(8) 

24(4) 

Ffoid 

11.7(4) 

13.3(5) 

12(3) 

Ffc 

6.69(2) 

6.70(3) 

6.68(5) 

O'. Fe 

0.28(3) 

0.27(4) 

0.24(7) 

a4 

Fe 

7.2(7) 

4.9(6) 

1.2(2) 

Ecyc 

4.3(2) 

4.3(1) 


°"cyc 

1.2(3) 

1.1(1) 


F ab 

10.1 

6.5 

1.2 

^unab 

13.0 

8.3 

1.6 

xL / dof 

0.94/549 

1.04/531 

1.06/487 


a in units 10 22 atoms cm -2 

b K a line normalization in units 10 _3 ph cm -2 s -1 

c Absorbed and unabsorbed fluxes in units 10 -9 erg cm -2 s -1 and in 0.1 - 
120 keV range. 


with high-energy cut-off (see Fig. [9]). We divided the data in 10 
equally spaced phase bins. We fixed the absorption component at 
the average value W H = 5.7^2x 10 22 atoms cm -2 . The emission iron 
line at ~ 6.7 keV and the absorption-like feature at ~ 4.5 keV were 
also included in the phase resolved analysis. We can see the clear 
anti-correlation of the photon index with phase flux. The parame¬ 
ters of the absorption-like and iron lines (line centroid and width) 
did not exhibit any significant phase variation. The cut-off and fold¬ 
ing energies appeared to be anti-correlated with each other. 


Table 2. Spectral parameters of the persistent spectrum of GRO J1744-28 
observed by BeppoSAX in 1997-03-21 for various spectral models. All en¬ 
ergies and line widths are given in keV. 


Par. 

PHC 

PHC 

no CRSF 

FDCO 

NPEX 

CompTT 

N a 

5.7(2) 

5.23(2) 

5.55(8) 

5.16(6) 

5.84(3) 

r 

1.26(7) 

1.147(3) 

1.09(4) 



Fcut 

18.4(1) 

18.1(1) 

15(1) 



Ffoid 

11.7(4) 

11.76(7) 

10.2(2) 

6.1(1) 


a\ 




0.48(4) 


CV2 




-2.0 


kT s 





0.18(3) 

kT e 





5.82(3) 

T e 





13.1(1) 

F Fc 

6.69(2) 

6.71(2) 

6.70(2) 

6.68(2) 

6.72(2) 

O'. Fe 

0.28(3) 

0.38(2) 

0.26(3) 

0.23(3) 

0.34(3) 

N b 

v Fe 

7.2(6) 

10.2(4) 

5.8(5) 

5.8(6) 

8.8(5) 

Ecyc 

4.3(2) 


4.47(9) 

4.55(5) 

4.2(3) 

C eye 

1.2(3) 


1.1(2) 

1.2(1) 

0.7(5) 

XrJ doi 

0.94/549 

1.28/552 

0.95/551 

1.04/551 

1.25/553 


a in units 10 22 atoms cm 2 

b K a line normalization in units 10 _3 ph cm -2 s -1 


3.2.1 Analysis of the burst spectra 

To investigate the burst spectra we aligned and stacked all observed 
bursts using the rising edge as a reference to improve statistics. 
We also subtracted the contribution of the persistent emission. We 
found that the shape of the combined spectrum of the bursts departs 
significantly from that of the non-bursting spectrum (Figj8]l. 

We first described the burst spectrum with the PHC and 
CompTT models. Both models give formally acceptable results (see 
Table|4j. The spectral param eters of the PHC model are similar to 
the ones reported recently bv lYounes et alJd2Q15h for the burst spec¬ 
tra (r = 0.2 + 0.1 and £ fo id = 7.6 + 0.5 keV). Such a hard power 
law together with a soft cut-off are rather peculiar and not typical 
of X-ray pulsars. This might suggest that the model describes a 
bump-like feature at soft energies. In addition, the best-fit results 
for the bursts spectrum using the CompTT model reveals a signif¬ 
icant change, with respect to the persistent spectrum, of both the 
seed photon temperature and the absorption column (seeTable 0]), 
which is difficult to understand. On the other hand, the apparent 
change of the seed photon temperature might suggest an additional 
soft component in the burst spectrum. Indeed, including a black- 
body component in either model allows to achieve comparable (or, 
in fact, slightly better) fit statistics while the other parameters of the 
continuum remain close to ones measured for the persistent spec¬ 
trum. In other words, we find that the burst spectrum significantly 
differs from the persistent one at soft energies, and this change is 
well accounted for by adding a soft blackbody-like component to 
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0.0 0.2 0.4 0.6 0.8 1.0 

Phase 


Figure 9. Spectral parameters of the pulsar GRO J1744-28 in 1997-03-21 
with the PHC continuum model as function of pulse phase. Phase: 0-0.2- 
0.4-0.6-0.8-1. Dashed lines are the MECS and PDS pulse profiles, they are 
the same. 


the unchanged continuum of the persistent spectrum. The black- 
body component has a temperature about 2.1 - 2.2 keY and a size 
close to a neutron star radius (for an assumed distance of 8 kpc). A 
Fe-K line is also marginally significant in the burst spectra of the 
observation on 1997-03-27. We have summarised the fit results for 
the 1997-03-21 data in Table 0 

To investigate the pulse phase dependence of the black body 
and power law components we carried out the phase resolved anal¬ 
ysis of the burst spectrum. Unfortunately, the statistics was insuffi¬ 
cient to constrain all model parameters, so we kept most of the pa¬ 
rameters fixed to the best-fit phase average values, and only allowed 
the normalizations of the continuum components to vary. The re¬ 
sults are shown in Fig. |T0| The pulse fraction for the soft and hard 
components are 6.5% + 3.4% and 13.8% + 1.9% respectively, i.e. 
the hard component varies with pulse-phase stronger than the soft 
one. 

To explore the time evolution of both components along the 
burst, we performed also time-resolved spectral analysis using the 
stacked data of all bursts. It is interesting to note, that both the tem¬ 
perature and normalization of the blackbody component (and thus 
the size of the emitting region) change along the burst as shown in 
Figure fill 



0.4 0.6 

Phase 


Figure 10. Changing normalizations of the soft (top) and hard (bottom) part 
of the bursts spectrum with phase for the GRO J1744-28 in 1997-09-21. 
lies in region from 0.97 to 1.45. 




Figure 11. Cooling of the soft component, during the bursts. The probability 
that the observed temperatures are due to a statistical fluctuation is < 1.7 x 
10 -8 (from Kolmogorov-Smirnov test), so the trend is significant. 
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Table 4. Best-fit parameters of the bursts spectrum GRO J1744-28 in 1997- 
03-21. All the energies and the line widths are given in keV. 


Parameter 

PHC 

bb+PHC 

CompTT 

bb+CompTT 

N a 

5.0(3) 

5.1(7) 

3.4(2) 

5.2 

kT bb 


2.1(3) 


2.2(1) 

R bb 


~ 7.5 


~ 7.6 

r 

0.26(8) 

0.8(4) 



Ecut 

4(1) 

18(5) 



Efold 

7.6(2) 

9(1) 



kT s 



1.51(8) 

0.2 

kT e 



5.4(2) 

5.2(1) 

T e 



4.6(5) 

18(1) 

*L/ dof 

0.96/408 

0.95/406 

1.03/409 

1.02/409 


a \n units 10 22 atoms cm -2 , ^in km. The absorbed and unabsorbed fluxes in 
0.1 - 120 keV energy range are F a b = 6.09 x 10 -8 erg cm -2 s -1 , F unab = 
6.83 x 10 -8 erg cm -2 s -1 . 


4 DISCUSSION 


The non-bursting spectrum of GRO J1744-28 has been measured 
for the first time in a broad energy range 0.7 to 120 keV with Bep- 
poSAX during the source outburst in 1996-1997. The persistent 
spectrum of the source is similar to the one of accreting pulsars 
and can be described by phenomenological models typically used 
for pulsars, e.g. by an absorbed power law with cut-off at around 
~ 18 keY. In addition, a line-like absorption feature at E ~ 4.5 keV 
is observed and required regardless of the used continuum model. 
We incidentally observe that there are no systematic effects known 
for the MECS and LECS, in the energy range where residuals are 
observed. In addition, this absorption-like feature is not observed in 
the spectrum of the bursts, which has similar shape and comparable 
statistical quality. Therefore, the feature is unlikely to be an artefact 
of the continuum modelling. 

Cyclotron resonance scattering features (CRSFs), interpreted 
as a result of resonant scattering of continuum photons off plasma 
electrons in strong magnetic fields close to the electron gyro fre¬ 
quency, are often observed in the spectra of accreting pulsars. 
CRSFs appears as a flux suppression of the continuum at cyclotron 
energies related to the magnetic field strength in the line forming 
region: 


E cyc * 11.57 B \2 (1 + z) -1 keV, 


where B 12 is the magnetic field strength in the units of 1 0 12 G , 
z is the gravita t ional r edshift (see Pottschmidt et al] (|2012|) . 
iDoroshenko et al.l d201C)h . ISantangelo et al.l ( 19991) for a recent re¬ 
view). By interpreting an absorption-like feature observed at E ~ 
4.5 keV in the spectrum of GRO J1744-28 as a CRSF, a magnetic 
field of B ~ 3.7 x 10 11 (1 + z) G is obtained for the neutron star of 
the binary system. This is in strikingly good agreement with earlier 
estimates of the magnetic field of this pulsar, which fall in range 
2-7 xlO 11 G (Finger et all 19961: Rappanort & Jossl[l997l : ICuill 1997b 


and (Younesetal 


In particular, 


20151; ID Ai et al 


D’Ai et al] d2015h have recently reported the de¬ 


tection of such an absorption feature at ~ 4.7 keV based on XMM- 
Newton and INTEGRAL data, which they also interpret as a cy¬ 
clotron resonance scattering feature (CRSF). Our result is in ex¬ 
cellent agreement with these findings. 

Besides the lack absence of absorption feature, the burst spec¬ 
trum was found to differ significantly from the non-bursting spec¬ 


trum, particularly at soft energies. The difference is most easily ac¬ 
counted for by inclusion of a black-body component on top of the 
power law-like persistent spectrum of accreting pulsars. The flux 
in 0.7 - 120 keY energy range measured for the power law and 

s" 

-2 c -1 


the blackbody components is F bb « 1.6 X 10 8 erg cm 2 s 1 and 


Fpower ~ 4 x 10“ 8 erg cm -2 s _1 , respectively. This indicates that 
~ 30% of the energy released in bursts arises from the soft thermal 
component not present in the persistent emission. 

The nature of the soft component is unclear and several hy¬ 
potheses may be considered. First of all, spectra of accreting pul¬ 
sars are known to change with luminosity and the luminosity in 
GRO J1744-28 increases by an order of magnitude during the 
bursts, so some intrins ic variation of the spectrum might be antic¬ 
ipated. For instance, iReig & Nespolil d2013h reported evident spec¬ 
tral softening for sever al pulsars above the so-called critical lu¬ 
mino sity ~ 10 37 ergs _1 (Bask o & Sunvaevlll976l : iMushtukov et al.l 
120151) . and similar softening might occur in GRO J1744-28, which 
at L x ~ 10 38 erg s -1 is likely also above this limit. 

Emission from the accretion column expected to form at this 
stage can also irradiate the accre tion dis k or the surface_ of neu¬ 
tron star (iFvubarskii & Svunvaevlll988l : IPoutanen et al.ll2013fo . and 
the thermalization of intercepted emission could be responsible for 
the soft component as well. The effective temperature of the irra¬ 
diated surface r eff « (L x /4^t 2 /ct S b) 1/4 ~ 2 keY is, in fact, com¬ 
parable with the observed. On the other hand, strong gravitational 
separation of ch emical elements in the neutron star atmospheres 
(lHameurv et al.l Il983h implies that the upper atmospheric layers 
consist almost exclusively of fully ionized hydrogen plasma which 
reflects most of the incident flux due to electron scattering. There¬ 
fore, the fraction of thermalized em ission is like ly to be less than 
10% (see a detail argumentation in IPoutanen et al.112013b . in con¬ 
trast with the ~ 30% observed. It is also hard to explain the ob¬ 
served pulsation of the soft component if it originates in the inner 
regions of accretion disk, although the even softer com p onent with 
kT ~ 0.5 keY reported recently bv lYounes et al.l d2015b : lD’Ai et al.l 
d2015b might indeed be related to irradiation of the accretion disk 
by the pulsar. 

Finally, taking into consideration the observed cooling of the 
thermal component along the burst (see Fig . fTTI) very similar to that 
observed in classical bursters dFewin et all fl993l) thermonuclear 
flashes on the NS surface could be also responsible for the observed 
soft emission. In fact, several authors have already considered the 
possibility that some of the bursts in GRO J 1744^28 mightbe Type- 
I bursts. Based on the BATSE observations. iFewin et al.ldl99d ) con¬ 
cluded that this is unlikely because the amount of matter accreted 
between the bursts is insufficient to explain the observed burst flu- 
ence if bursts are powered by thermonuclear burning. Indeed, the 
inter-burst to burst fluence ratio a ~ 4 deduced from BATSE ob¬ 
servations was found to be much small er tha n a ^ 40 typical of 
thermonuclear hydrogen burning ( Fewin etldlll993l) . 

On the other hand, Ijahoda etldl (Il997h . based on the RXTE 
observations close to the peak of the outburst, found significanly 
higher a ~ 34. Moreover, we found that the flux of the thermal 
component constitutes only about 30% of the total flux and, and 
therefore of the burst fluence. If this is taken into account, even 
at lower luminosities such as during the BeppoSAX observations, 
when we estimate a ~ 5 - 15 for the bolometric flux, the same ratio 
calculated for the thermal component alone is factor of three higher, 
i.e. still comparable w ith that obse rved in classical bursters. There¬ 
fore the argument by IFewin et all ( 1996 ) clearly does not always 
hold. 

The strongest argument against the thermonuclear origin of 
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the soft component in burst spectrum of GRO J1744-28 is that 
it is an accreting pulsar. Indeed, the temperature and pressure 
of matter funneled by the magnetic field to the polar areas are 
usually sufficient for stable thermonuclear burning, so no un¬ 
burned matter accumulates at the surface and thus no Typ e-I bursts 
are ob served from accreting pulsars. In fact, iBildsten & Brownl 
dl997l) did consider the possibility of unstable thermonuclear 
burning in GRO J1744-28, and concluded that under the most 
favourable assumptions it should not be possible. In particular, 
iBildsten & Brownl ( 1997 ) argue that to escape the accretion flow 
and spread over the NS surface (where it can be accumulated and 
subsequently ignite as a thermonuclear flash), the accreted plasma 
must overcome the magnetic pressure and turn the magnetic field 
lines parallel to the neutron star surface. This implies a pressure of 
P ~ 10 24 erg cm -3 while the hydrogen/helium mixture burns in the 
stable regime at P « 10 22 erg cm -3 . 

In other words, accreting matter remains confined at the polar 
caps where it burns ste adily as expected. Howeve r, in our opinion, 
this key assumption of IBildsten & Brown (1997 ) is probably too 
conservative, and, in fact, corresponds to the case when the plasma 
is not confined to polar areas at all. This would basically imply that 
GRO J1744-28 is an ordinary burster, which clearly it is not. On 
the other hand, the magnetic pressure at the poles of the neutron star 
can be estimated based on the observed CRSF energy and turns out 
to be B 2 /8n « 10 22 erg cm" 3 comparable with the pressure required 
for stable thermonuclear burning. Therefore, the magnetic field is 
unlikely to confine plasma at significantly higher pressures and part 
of accreted matter probably leaves the polar areas before burning, 
thereby creating conditions for thermonuclear flashes. We note also 
that, if realised, the onset of a such flash could in principle disrupt 
the inner parts of the accretion disk triggering the instabilities re¬ 
sponsible for the enhanced accretion rate and most of the observed 
burst emission. 

In conclusion, we fully agree that there is no doubt that the 
mass accretion rate does increase during the bursts GRO J1744-28 
and is responsible for the bulk of the observed emission. However, 
some fraction of the flux in principle could still be due to unstable 
thermonuclear burning at the NS surface. 


5 SUMMARY AND CONCLUSIONS 

We presented the results of the analysis of three BeppoSAX obser¬ 
vations, with a total exposure time 270 ks, carried out in the de¬ 
clining phase of the 1997 outburst of the unique bursting pulsar 
GRO J1744-28. 

Pulsations with a period of 0.4670 s with a stable pulse profile 
were detected in all observations. The pulsed fraction was found 
to vary with energy, reaching a minimum of 10-20% in the en¬ 
ergy range 5-40 keV and increasing for higher and lower energies, 
especially at lower luminosities. Several tens of bursts with typi¬ 
cal durations of about 10 s were detected as well. Depletion in the 
light curve is observed after at least some of the bursts. The source 
luminosity typically increases by factor of ten during the bursts, al¬ 
though a number of dimmer bursts are also observed. 

The non-bursting broad band X-ray spectrum was found to be 
well described by several phenomenological models typically used 
for accreting pulsars. An iron line at ~ 6.7 keV and an absorption 
feature at ~ 4.5 keV were also required to fit the data. We inter¬ 
pret the absorption feature as a cyclotron line which implies a mag¬ 
netic field of the neutron star of B ~ 3.7 x 10 11 G. This value is 
in good agreement with earlier predictions and with the recent ob¬ 


servations of the 2014 outburst with XMM-Newton and INTEGRAL. 
iD’Ai et all (120151) have detected the same feature at E ~ 4.7 keV 
in agreement, within uncertainties, with our measurement. 

The average burst spectrum could be represented as a combina¬ 
tion of the harder non-bursting spectrum and an additional soft ther¬ 
mal component with temperature of about 2 keV. The burst spec¬ 
trum requires that neither the iron 6.7 keY line nor the absorption 
feature at 4.5 keV. Both components are pulsed although the ampli¬ 
tude is smaller for the soft component, ~6.5%vs~13.8%. 

We discussed a possible nature of the thermal component and 
speculate that it could be caused by thermonuclear flashes which 
possibly trigger the accretion rate enhancement responsible for the 
bursts observed in the source. This hypothesis is based on several 
similarities between typical Type I bursts and the thermal compo¬ 
nent, namely, the burst duration of 10 s, the inter-burst to burst flu- 
ence ratio and cooling as a burst progresses. There are some ob¬ 
servational and theoretical arguments that disfavor this hypothesis. 
These cannot be verified with the existing data, however, addit ional 
data from current or future missions like LOFT din ’t Zand et al.l 
2015) might provide decisive insights on this puzzling source. 
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